A fractional-order mathematical model for COVID-19 outbreak with the effect of symptomatic and asymptomatic transmissions

The purpose of this paper is to investigate the transmission dynamics of a fractional-order mathematical model of COVID-19 including susceptible (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsc {S}$$\end{document}S), exposed (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsc {E}$$\end{document}E), asymptomatic infected (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsc {I}_1$$\end{document}I1), symptomatic infected (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsc {I}_2$$\end{document}I2), and recovered (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsc {R}$$\end{document}R) classes named \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {SEI_{1}I_{2}R}$$\end{document}SEI1I2R model, using the Caputo fractional derivative. Here, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {SEI_{1}I_{2}R}$$\end{document}SEI1I2R model describes the effect of asymptomatic and symptomatic transmissions on coronavirus disease outbreak. The existence and uniqueness of the solution are studied with the help of Schaefer- and Banach-type fixed point theorems. Sensitivity analysis of the model in terms of the variance of each parameter is examined, and the basic reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(R_{0})$$\end{document}(R0) to discuss the local stability at two equilibrium points is proposed. Using the Routh–Hurwitz criterion of stability, it is found that the disease-free equilibrium will be stable for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{0} < 1$$\end{document}R0<1 whereas the endemic equilibrium becomes stable for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_{0} > 1$$\end{document}R0>1 and unstable otherwise. Moreover, the numerical simulations for various values of fractional-order are carried out with the help of the fractional Euler method. The numerical results show that asymptomatic transmission has a lower impact on the disease outbreak rather than symptomatic transmission. Finally, the simulated graph of total infected population by proposed model here is compared with the real data of second-wave infected population of COVID-19 outbreak in India.


Introduction
Mathematical models are powerful tools for understanding and preventing the infectious disease transmission. The researchers gain valuable information for several contagious diseases by studying stochastic and deterministic models. Kermack and McKendrick (1927) created a helpful model for executing and growing complex epidemic models, now regarded as a fundamental model in an epidemiology study [1]. Many infectious diseases can transmit in both vertical and horizontal paths. HIV/AIDS, Herpes Simplex, Hepatitis B, and Rubella are a few examples of such infections. These diseases are horizontally transmitted in people and animals by host-to-host interaction or disease carriers such as flies, mosquitoes, and others. The reader can find several good papers on the infectious disease modeling in the literature, for example, HIV/AIDS [2], Zika [3], Ebola [4], West Nile [5], Influenza [6], Dengue [7], Oncolytic [8], Cancer [9].
The first epidemic of severe acute respiratory syndrome (SARS) occurred in mainland China in 2003, while another outbreak is known as Middle East respiratory syndrome (MERS) occurred in South Korea in 2015. The new virus (2019-nCoV), which is highly transmissible and virulent, was identified in a single individual in the Chinese city of Wuhan [10]. This new disease, which produces a severe acute respiratory syndrome, has spread throughout the world. There have been over 253 million confirmed COVID-19 registered cases, and over 5.1 million fatalities globally since November 10, 2021 [11]. America, Europe, Africa, Southeast Asia, the Western Pacific, the Eastern Pacific, and the Mediterranean are the most effected by the coronavirus. The first signs of a COVID-19 infection include a dry cough, fever, exhaustion, and shortness of breath, which develop in 2-10 days and can lead to pneumonia, SARS, kidney failure, and even death [12]. The pandemic continued to expand, but vaccinations have now halted the transmission. Many researchers have been developed mathematical models for the transmission of the aforesaid disease involves various compartments, for example Susceptible-Infected (SI) [13], SI-Removed (SIR) [14], S-Exposed-IR (SEIR) [15], SEI-Quarantined-R (SEIQR) [16], and SIQ-Home Quarantined-R-Diseased (SIQHRD) [17].
For a year, there were no particular treatment approaches such as antibiotics or other treatments against the coronavirus because it is a single RNA virus. As a result, non-pharmaceutical measures such as isolation, social distance, quarantine, awareness programs, developed an SEI 1 I 2 R model for COVID-19 disease transmission, in which the population is divided into five classes: susceptible S, exposed E, asymptomatic infected I 1 , symptomatic infected I 2 , and recovered or removed R. The model assumes the total population N is not constant throughout the time and equal to the sum of all compartment sizes at any time t, i.e., N = S + E + I 1 + I 2 + R. We made the following assumptions to formulate the model: Asymptomatic transmission is represented by the primary pathway SEI 1 R, in which asymptomatic infected individuals spread the disease at a rate of μ 1 . The SEI 2 R is the second pathway that denotes indicative transmission, in which symptomatic individuals transmit the infection at a rate of μ 2 . The symptomatic case's transmission rate μ 2 and infection rates ρ 2 are considered to be higher than the asymptomatic case's, i.e., μ 2 > μ 1 and ρ 2 > ρ 1 . Asymptomatic patients have a higher recovery rate than symptomatic patients (β 1 > β 2 ), and asymptomatic patients have a lower mortality rate (α 2 > α 1 ) than symptomatic patients. is the recruitment rate, representing the increment of compartment S in the form of birth or migration or recovered individuals. The natural death rate α is equal to the recruiting rate into compartment S. In this study death population due to the COVID-19 infection is only counted.
where μ i , ρ i , β i , and α i (i = 1, 2) denotes the transmission rate, infection rate, recovery rate, and death rate in the mentioned above pathways, respectively. λ 1 and λ 2 are the positive constants. The quantities μ 1 S I 1 N and μ 2 S I 2 N in the first equation represent the interaction of susceptible individuals with asymptomatic and symptomatic infectious peoples with μ 1 and μ 2 rates, respectively. As can be seen in the second equation, these interactions increase the number of people who are exposed. The terms ρ 1 E and ρ 2 E denote the number of people becoming infected via the exposure of susceptible people. In the third equation, α 1 I 1 and β 1 I 1 show the number of death and recovered population in the asymptomatic class. The terms α 2 I 2 and β 2 I 2 in the fourth equation specify the number of death and recovered individuals in symptomatic class. The term αS represents the natural death of susceptible people and αR shows the death of recovered people. Dividing all the equations of the model (1) by (N) and setting S = S N , E = E N , , as a result, model (1) may be expressed as: where S + E + I 1 + I 2 + R = 1. The proposed mathematical model for COVID-19 outbreak is considered under the initial conditions: nonnegative. The model is designed to look at the role of asymptomatic transmission rates in the COVID-19 epidemic in comparison to symptomatic transmission rates. Using the above assumptions, we propose the following SEIR compartmental model of order 0 < θ ≤ 1 at time t ∈ [0, T], T < ∞ as: subject to the following initial conditions S(0), E(0), I 1 (0), I 2 (0), R(0) ≥ 0. The flowchart for the proposed model (3) can be seen in Fig. 1.

Preliminary definitions
This section is dedicated to some basic definitions and lemma that are necessary for achieving the existence theory.
Definition 3.1 [27] The integral of order θ > 0 of a function Z ∈ L 1 ([0, T], R + ) in the sense of R-L is the following: Definition 3.2 [27] For derivative of order θ > 0 of a function Z ∈ AC n (0, ∞), in the sense of Caputo is the following: Lemma 3.1 [27] If θ > 0, then the differential equation C D θ Z(t) = G(t) will have the following solution: where n = [θ ] + 1.

Existence theory of model (3)
We construct the existence theory for the considered model (3) in this section. Let us write model (3) as follows: The proposed problem (3) can be reformulated in the following form: S, E, I 1 , I 2 , R), S, E, I 1 , I 2 , R), S, E, I 1 , I 2 , R), The considered system (3) can be written in the form presented using Eq. (4): Equation (6) is solved using the R-L type integral as follows: where For the existence theory, we define Banach space under the norm Z = (S, E, I 1 , Let U : J → J be an operator defined as follows: The existence of a solution is proved using the following theorem.  (P) There is constant L G > 0 such that for Z, Z ∈ J, then is hold. There is at least one solution for the considered system (3) .
Proof We begin by proving that the operator U is completely continuous. Suppose a sequence {Z n } such that Z n → Z in J, then for t ∈ [0, T], we have Then by definition of S, Then for each Z ∈ S, we can obtain Thus, U is uniformly bounded.
Further, suppose 0 ≤ t 2 ≤ t 1 ≤ T. Then As a result, U is equicontinuous. Since U is both continuous and bounded, it is compact and so completely continuous. Let = Z ∈ J : Z = U(Z), ∈ [0, 1] , we need to confirm that is bounded. Suppose Z ∈ , then for t ∈ [0, T], we have: Hence, the operator is completely continuous and the set is bounded. Therefore, according to Theorem 4.1, U has at least one fixed point. As a consequence, the considered system (3) has the same number of solution.
To prove the uniqueness, we utilize the Banach's fixed point theorem [56].

Theorem 4.3 If the assumption (P) is hold and
Now, in view of (P), for Z, Z ∈ J, and for each t ∈ [0, T], we have The contraction of U proves that the solution to the problem under consideration is unique.

Well-posedness and biological feasibility
This section looks at the interval and region where our system's solution will make perfect historical sense. Let N be the net population size, i.e., Then Since, R ≤ N. So the above inequality becomes: We want that the function N(t) is a positively increasing function N (t) > 0 then N(t) < α and N (t) ≥ N(t). The aforementioned inequality is known as the threshold population level in the literature. As a result, we infer that the accepted set of solutions to the proposed model must be contained within: Here in biological terms R 5 + is the positive cone of R 5 that also contains its lower-dimensional faces.
6 Disease-free and endemic equilibrium point The point that no disease exists in the population is known as the disease-free equilibrium point.
Applying the required conditions E = I 1 = I 2 = 0 in the system (12) for disease-free equilibrium, then the disease-free equilibrium point of the considered model (3) is given below:

Basic reproduction number
The virus's reproduction number, also known as R 0 , is a way of rating a disease's propensity to spread and calculate the average number of individuals infected by one infected person. In general, this basic reproduction number R 0 is an epidemiologic metric used to describe the transmissibility of infectious agents or contagiousness. So in this section, we will utilize the next generation matrix technique to calculate the basic reproductive number for this model. Therefore, the infected compartments of our model are the following: The Jacobian matrix of the system (14) is given by: Now decomposing the matrix J 0 in terms of A and B, that is J 0 = A − B, we have: The basic reproduction number is obtained from R 0 = ρ(AB −1 ), so we have: where S 0 = α . Hence, Eq. (15) is the basic reproduction number for the proposed system.

Stability of equilibrium points
The stability of equilibrium points is investigated in this section. The Jacobian matrix of system (3) is as follows: So, the Jacobian matrix of system at E 0 is: Theorem 8.1 Let R 0 and A i , (i = 1, 2, 3) be defined by (15) and (17), respectively. If R 0 < 1, then by Routh-Hurwitz stability criterion [57], the disease-free equilibrium point E 0 of the model (3) is locally asymptotically stable if and only if A 1 > 0, A 3 > 0 and A 1 A 2 > A 3 .
Proof The characteristic equation of the Jacobian matrix at the disease-free equilibrium point is det(J(E 0 ) − φ I ) = 0. Then one can achieved the following: where ⎧ ⎪ ⎨ ⎪ ⎩ Rewriting the A 3 in terms of the basic reproduction number R 0 , A 3 becomes: It is obvious in the first factor of above characteristic Eq. (16) that two eigenvalues φ 1 and φ 2 are always negative for α > 0. The other three eigenvalues can be obtained by solving the following equation: The remaining three eigenvalues have negative real parts, if they follow the Routh-Hurwitz criteria of order 3, such that A 1 > 0, A 3 > 0 and A 1 A 2 > A 3 . Therefore, R 0 < 1 if and only if all the eigenvalues have negative real parts. Which completes the proof.
Individuals who are asymptomatic yet infectious play a critical role in spreading the disease across the population. When R 0 > 1, the endemic equilibrium state of the model (3) is stable, as shown in the following theorem.
Then by Routh-Hurwitz stability criterion [57], the endemic equilibrium point E * of the model (3) (3) about E * such that: where The characteristic equation associated with J(E * ) is: where From Eq. (19) it is clear that there are five corresponding eigenvalues of J(E * ). One of them is φ 1 = −α, having negative real part. The remaining four can be obtained by the solution of the following equation: Using Routh-Hurwitz criterion of order 4, we can conclude that given system (3) is local asymptotic stable at E * .

R 0 Sensitivity analysis
The sensitivity analysis is applied to study the effect of the parameters on proposed COVID-19 outbreak model. In particular, it is necessary to identify the most sensitive parameters causing a disturbance in model dynamics by a small change in numeric value. To check the R 0 sensitivity, we calculate its derivatives as follows: Since, all the parameters are positive so It concludes that the reproduction number (R 0 ) is increasing with , μ 1 , μ 2 and ρ 2 . The normalized sensitivity indices corresponding to these parameters are estimated as follows: Here, the sensitivity index can be constant depending on some parameters or can be free of any independent parameters. The partial rank correlation coefficient (PRCC) results for significance of parameters involved in R 0 can be seen in Fig. 2 Table 2.

Numerical results and discussion
In this section, we discuss the impact of parameters on reproduction number R 0 in numerical simulations of the proposed model using the fractional Euler method [58]. The numerical simulation of the present model (3) will be carried out to support the analytical result using the value of the parameters as given in Table 1.
The sensitivity analysis is carried out by estimating the sensitivity indices based on the parameter values to examine the parameter impact on R 0 . The sensitivity indices given in Table 2 are obtained by using the values of involved parameters. From the computed sensitivity indices, it can be seen that a 10% increase in the recruitment rate , asymptomatic transmission rate μ 1 , symptomatic transmission rate μ 2 , and symptomatic infection rate ρ 2 cause to increase the value of R 0 by 10%, 2.6%, 7.3%, and 1.8%, respectively, and can lead to an outbreak subsequently. However, it can be noticed that the symptomatic infection rate ρ 2 does not affect R 0 significantly. But on the other hand, asymptomatic infection rate ρ 1 , asymptomatic recovery rate β 1 , symptomatic recovery rate β 2 , asymptomatic death rate α 1 and symptomatic death rate α 2 describe that increasing their values by 10% will decrease the value of R 0 by 1.8%, 2.5%, 5.5%, 0.1%, and 1.7%, respectively.
The impact of some parameters on the reproduction number is shown in Fig. 3. Figure 3a describes the reproduction number R 0 sensitivity versus the asymptomatic transmission rate μ 1 and asymptomatic death rate α 1 . Since from Table 2, it is clear that μ 1 has a positive impact while α 1 has a negative impact on reproduction number R 0 . So, it can be seen in Fig. 3a the asymptomatic death rate has no impact on R 0 because no variation can be seen in R 0 when α 1 increases. On another side, the asymptomatic transmission rate μ 1 has a very high impact on R 0 . Figure 3b describes the reproduction number R 0 sensitivity versus symptomatic infection rate   −0.55385 ρ 2 and asymptomatic transmission rate μ 1 . So it can be seen in Table 2 that both parameters has a positive impact on reproduction number R 0 . But the impact of asymptomatic transmission rate is higher than the symptomatic infection rate, which is depicted in Fig. 3b. Figure 3c shows the reproduction number R 0 sensitivity versus the symptomatic transmission rate μ 2 and asymptomatic death rate α 1 . The impact of the symptomatic transmission rate is positive on the reproduction number, while the impact of the asymptomatic death rate is negative. The impact of asymptomatic death rate is almost negligible, but the effect of the symptomatic transmission rate is quite high enough where can be seen in Fig. 3c. Figure 3d shows the reproduction number R 0 sensitivity versus asymptomatic (a) (b) (c) (d) Fig. 4 The solution behavior of the all compartments for various values of θ recovery rate β 1 and symptomatic recovery rate β 2 . Both parameters have a negative impact on the reproduction number R 0 (see Table 2), but the affect is almost the same, which is shown in Fig. 3d. Figure 3e demonstrates the reproduction number R 0 sensitivity versus the asymptomatic transmission rate μ 1 and symptomatic recovery rate β 2 . The transmission rate μ 1 of the asymptomatic pathway has a positive impact on reproduction number R 0 and symptomatic recovery rate β 2 has a negative impact on R 0 , which can be seen in Table 2. From Fig. 3e, we can conclude that the symptomatic transmission rate μ 1 effects the reproduction number more than the symptomatic recovery rate β 2 . Figure 3f describes the reproduction number R 0 sensitivity versus symptomatic death rate α 2 and symptomatic recovery rate β 2 . Both parameters have a negative impact on reproduction number R 0 . But the impact is almost the same, which is shown in Fig. 3f. Now we discuss the numerical results of the governing model in terms of approximate solutions. To achieve this goal, we used the efficient fractional-order Euler technique with the Caputo fractional operator. The initial conditions are considered as S(0) = 1 − E(0) − I 1 (0) − I 2 (0) − R(0), E(0) = 5/20000, I 1 (0) = 0, I 2 (0) = 0, R(0) = 0 with the total population N = 1 and the values of the parameter are given in Table 1. We have utilized the MATLAB software to simulate the numerical results. Graphs are depicted the approximate solutions obtained by using the suggested numerical approaches under the value of the considered parameters. The dynamics of the suggested model population are graphically shown in Fig. 4 using the Caputo fractional derivative with time sequence framework. Figure 4a shows the population behavior for θ = 0.7, Fig. 4b shows the population behavior for θ = 0.8, Fig. 4c shows the population behavior for θ = 0.9 and the population behavior for θ = 1.0 is shown in Fig. 4d. Figure 5 shows the solution behavior of S(t), E(t), I 1 (t), I 2 (t), and R(t) for θ = 0.7, 0.8, 0.9, 1.0. Figure 5a shows that the susceptible population decreases rapidly at different rates and stabilizes after about 100 days due to restriction of individual movements. As we can see that the decay rate is faster for a small value of θ , since the susceptible population S(t) is decreasing, which affects the exposed population E(t) and more individuals get the infection. That is why we see an increase in the exposed population in Fig. 5b. Besides, the exposed population E(t) in Fig. 5b is increased very quickly after 15 days and a decrease occurs Combined effect of μ 1 and μ 2 , Combined effect Moreover, the exposed population increases and decreases rapidly for smaller values of θ as compared to bigger values. The asymptomatic class I 1 (t) and symptomatic class I 2 (t) in Fig. 5c and d have almost the same pattern, but the asymptomatic population becomes stable after 150 days and the symptomatic population becomes stable after 200 days. Both populations have grown rapidly, but the decline occurred due to standard operating procedures (SOPs) and vaccinations. The infection increases and decreases quickly for smaller values of θ . Since we know that many innocent people have been died due to this pandemic, but there are many who have survived as well due to timely treatment and following the SOPs. Most people have been vaccinated and are now protected from infection. 52.9% population have received their first dose of the COVID-19 vaccination worldwide till November 13, 2021 [61], which reduced the infection and that is why the recovery rate has almost been stable. So in Fig. 5e, it can be seen that the recovery rate is almost 83% of the population. The increase in recovery is due to a strong protection rate, people's awareness about the infection, government action, and vaccination. Besides, this rise up of recovery is also different for various values of fractional order.
Here, we discuss the impact of the parameters , μ 1 , μ 2 and ρ 2 on the basic reproduction number R 0 , where the values of mentioned parameters is given in Table 2. The sensitive indices (mentioned above parameters) will be used to assess the propensity of the total infection (I 1 + I 2 ) curve. Table 3 shows the loop iterations used to simulate model (3). The simulation has used different sensitive parameter values to generate graphs of total infection and basic reproduction number R 0 to examine the COVID-19 epidemic pattern while the other parameters remain unchanged. The graphs are displayed via a panel in Figs. 6, 7, 8, 9, 10, 11, where the left showing the basic reproduction number pattern and the right indicating the total infection behavior concerning a specific sensitive parameter.
According to the sensitivity analysis, the symptomatic transmission rate μ 2 is the second most sensitive component impacting the epidemic. The parameters μ 1 and μ 2 have nearly a similar effect on basic reproduction number R 0 . In comparison with asymptomatic transmission rate μ 1 in Fig. 6a, with symptomatic transmission rate μ 2 in Fig. 7a, μ 2 causes a quicker and higher increase in the total infection. The graph in Fig. 7b shows that roughly 58 percent of the overall population becomes infected with the In contrast, about 52 percent of the total population becomes infected with asymptomatic transmission rate μ 1 on day 48 in Fig. 6b. Consequently, maintaining self-isolation and quarantine would be the most effective ways to control COVID-19, which is disseminated by symptomatic transmission. Figure 8 depicts the effect of the symptomatic infection rate ρ 2 on basic reproduction number R 0 and on total infection. From Table 3 and Fig. 8a, it can be seen that the symptomatic infection rate ρ 2 is less sensitive to the R 0 . Moreover, as shown in Fig. 8b, it has a significant impact on the total infection of COVID-19 pandemic. Figure 9b shows the high effect of recruitment rate on R 0 . Also, it is shown for a minimal increase in , the total number of infected people increases significantly.
Imposing a ban on population migrating from one place to another would thus be the most effective way of handling this outbreak. The initial strategic action policy should be to established a management plan to control in order to keep R 0 < 1. As it is shown in Fig. 9a, the value of R 0 reaches from 5.8 to 20 for a minimal change in from 0.00004 to 0.00016, indicating that the recruitment rate to a region during the transmission of COVID-19 has a significant impact on the pandemic. Furthermore, Figs. 10 and 11 present the simultaneous impact of ( , μ i ), (i = 1, 2) on basic reproduction number and total infected population. Figure 10 illustrates the impact of ( , μ 1 ) on R 0 and the total infected population while Fig. 11 demonstrates the impact of ( , μ 2 ). Both total infection graphs 10b and 11b are almost identical but ( , μ 2 ) significantly impacts basic reproduction number compared to ( , μ 1 ). This means that the severe endemic will happen if the recruitment rate and symptomatic transmission rate change simultaneously.
Recently, India has faced the second wave of the COVID-19 pandemic, which had serious implications in the form of spiraling cases, restricted availability of essential medications, and increased mortality, particularly among the young population. After the second wave, the country was infected approximately by 34  November 13, 2021. The second wave was spread rapidly and had severely impacted the whole country. It is thought to be the second worst epidemic to hit India nearly after a century. The last menace occurred during the 1918 influenza outbreak, which claimed the lives of millions people. Here, we have taken the real data of COVID-19 spread in India from February 22, 2021, until November 10, 2021, from reference [62] and compared it with the simulated data in this research, which is shown in Fig. 12.

Conclusion
This study was aimed to investigate the effect of symptomatic and asymptomatic transmissions on the coronavirus (COVID-19) outbreak using a fractional-order mathematical model, which is achieved successfully. The existence theory and uniqueness of proposed model solution have been carried out with the help of Schaefer-and Banach-type fixed point theorems. The local stability of the model is established through the Routh-Hurwitz criteria, which is concluded that the disease-free equilibrium is stable for R 0 < 1, whereas the endemic equilibrium becomes stable for R 0 > 1 and unstable otherwise. The effect of the various parameters on the reproduction number (R 0 ) has also been analyzed. Besides, the transmission of the disease under various fractional orders is investigated. In Fig. 11, it can be seen that the recruitment rate and symptomatic transmission rate are significantly impacting the reproduction number. Therefore, the COVID-19 epidemic can be significantly reduced by controlling people from migrating and strictly enforcing personal measures simultaneously. The total infected population plot has been compared with the real infected population of India's second wave of COVID-19 outbreak (proportion of the population) in Fig. 12 which showed the great agreement of simulated results in this research with the real data. It is worth mentioning that the proposed model here can be studied under stochastic or fuzzy differential operators in future research.